function out = frameDiv(y)
[M N P] = size(y);
N=N/2;
p1 = y(1:M,1:N,:);
p2 = y(1:M,N+1:2*N,:);

out1(1:M,2:N-1,:) = p1(1:M,2:N-1,:)-p1(1:M,1:N-2,:);
out1(1:M,N,:) = -p1(1:M,N-1,:);
out1(1:M,1,:) = p1(1:M,1,:);

out2(2:M-1,1:N,:) = p2(2:M-1,1:N,:) - p2(1:M-2,1:N,:);
out2(M,1:N,:) = -p2(M-1,1:N,:);
out2(1,1:N,:) = p2(1,1:N,:);
out = out1 + out2;
end


